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ABSTRACT 

Context. This is the fourth in a series of papers that deal with angular momentum transport by internal gravity waves in stellar interiors. 
Aims. Here, we want to examine the potential role of waves in other evolutionary phases than the main sequence. 
Methods. We study the evolution of a 3 M Population I model from the pre-main sequence to the early- AGB phase and examine 
whether waves can lead to angular momentum redistribution and/or element diffusion at the external convection zone boundary. 
Results. We find that, although waves produced by the surface convection zone can be ignored safely for such a star during the main 
sequence, it is not the case for later evolutionary stages. In particular, angular momentum transport by internal waves could be quite 
important at the end of the sub-giant branch and during the early-AGB phase. Wave-induced mixing of chemicals is expected during 
the early-AGB phase. 
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1. Introduction 

In recent years, several authors studied the impact of internal 
gravity waves (IGWs) in a variety of main sequence stars. These 
waves were initially invoked as a source of mixing in stellar inte- 
riors in low-mass stars with an extended surface convection zone 
(Press 1981; Garcia Lopez & Spruit 1991; Schatzman 1993; 
Montalban 1994) and also as an efficient process in the synchro- 
nization of massive binary stars (Goldreich & Nicholson 1989). 
More recently, it was suggested that IGWs may play a role 
in braking the solar core (Schatzman 1993; Zahn, Talon, & 
Matias 1997; Kumar & Quataert 1997). This idea was con- 
firmed first in static models (Talon, Kumar, & Zahn 2002) and 
recently in the complete evolution of solar-mass models, evolved 
all the way from the pre-main sequence to 4.6 Gy (Charbonnel 
& Talon 2005). 

All these authors find that IGWs are easily excited, and a 
similar conclusion is reached by studies of convection on top of 
a stably stratified layer in 2-D and 3-D hydrodynamic numeri- 
cal simulations (e.g.Hurlburt et al. 1986, 1994; Andersen 1994; 
Nordlund et al. 1996; Kiraga et al. 2003; Dintrans et al. 2005; 
Rogers & Glatzmeier 2005a,b). These waves have a strong im- 
pact on stellar evolution, especially through their effect on the 
rotation profile. Through differential filtering, IGWs indeed play 
a major role in the redistribution of angular momentum in stars, 
which determines the extent and magnitude of rotation-induced 
mixing (Talon & Charbonnel 2005, hereafter TC05). 

In the spirit of applying a unique set of physical principles to 
stellar evolution, one needs to assess the impact of such waves 



Send offprint requests to: Suzanne Talon 



on stars of all masses and at various evolutionary stages. It is our 
purpose in this series of papers to consistently examine the full 
Hertzsprung-Russel diagram (HRD) for determining when such 
waves are efficiently emitted and how they can affect stars when 
their complete rotational history is being considered. 

In Talon & Charbonnel (2003, hereafter Paper I; see also 
Talon & Charbonnel 1998), we showed how the appearance of 
IGWs in solar-metallicity main sequence stars with an effec- 
tive temperature T e ff < 6700 K (i.e., when the surface convec- 
tion zone becomes substantial) can explain the existence of the 
lithium dipQ in stars undergoing rotational mixing. 

In Talon & Charbonnel (2004, hereafter Paper II), we exam- 
ined the IGW generation in Population II main sequence stars. 
We showed that, along the lithium plateatfl the net angular mo- 
mentum luminosity of IGWs is constant and high enough to en- 
force quasi solid-body rotation similar to that of the Sun in these 
stars. We proposed that this behavior could play a major role in 
explaining the constancy of the lithium abundance in the oldest 
dwarf stars of our Galaxy (see also a discussion in Charbonnel 
& Primas 2005). 

Here, we wish to look at other evolutionary stages, especially 
the pre-main sequence (PMS) and the more advanced stages for 
intermediate-mass stars. We focus in particular on stars in which 

1 The so-called lithium dip, observed in all galactic open clusters 
older than ss 200 Gyr, as well as in the field, refers to a strong lithium 
depletion in a narrow region of sa 300 K in effective temperature, cen- 
tered around T cff = 6700 K (e.g. Wallerstein et al. 1965; Boesgaard & 
Tripicco 1986; Balachandran 1995). 

2 The lithium plateau refers to the constant Li abundance of 
Population II dwarf stars above T cff ss 5500 K (e.g. Spite & Spite 1982). 
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Fig. 1. Evolution of our 3 M Pop I model in the HR diagram, 
from the PMS to the early-AGB. Numbers correspond to se- 
lected evolutionary points that are listed in Table Q] For clar- 
ity, the PMS is drawn with a dashed line, and the corresponding 
models are identified by squares, while the star symbol corre- 
sponds to the appearance of the convective core. The circles cor- 
respond to evolutionary points on the main sequence, sub-giant, 
and first ascent giant branches, while triangles refer to evolution- 
ary points on the clump and on the early-AGB. 

IGWs generation by the surface convection zone is limited on the 
main sequence, i.e., Pop I star originating from the left side of the 
Li dip. To do so, we follow the evolution of a 3 M , Z = 0.02 
star from the PMS up to the early asymptotic giant branch (early- 
AGB f| We estimate at which stages waves are efficiently gen- 
erated in the outer convection zone. This is a prerequisite for 
evaluating their impact on the rotation profile of the correspond- 
ing stellar models. We also look at the existence of a shear layer 
oscillation (or SLO) as a direct source of turbulence and mixing 
at the convective boundary. 

We begin in §|2]with a description of the evolution of relevant 
characteristics for our 3 M , Z = 0.02 star. In §[3] we recall the 
main aspects of the formalism we used to evaluate the impact of 
IGWs on stellar evolution. The following sections are devoted to 
discussing the pre-main sequence (§@), the main sequence and 
the sub-giant branch (§0, the red-giant branch (§|6]), and the 
clump and the early-AGB (§|7J- 

2. Evolution of relevant characteristics 

In this paper we focus on a 3 M Q with a metallicity Z = 0.02 
computed with the stellar evolution code STAREVOL (V2.30; 
Siess et al. 2000). Details regarding the code and the input 
physics may be found in TC05 and Palacios et al. (2006). The 
model was computed with classical assumptions, i.e., with nei- 
ther atomic diffusion nor rotation. Convection is treated using 
the standard mixing length theory with (Xmlt = 1-75. The treat- 
ment of convection has an impact on the excited wave spectrum 
(see § |3~TT >. 

When a star evolves, the characteristics that are relevant to 
wave excitation and to momentum extraction by IGWs change. 
We thus chose a series of models (i.e., evolutionary points) to 
correctly represent the evolution of the stellar structure in terms 

3 The peculiar case of thermally pulsing AGB (TP-AGB) stars will 
be presented elsewhere. 



Fig. 2. Evolution of the temperature at the boundaries of con- 
vection zones for the evolutionary points selected for the 3 M 
Pop I star and given in Table Q] Core and surface (i.e., T e ff) tem- 
peratures are also shown. The evolutionary phases are indicated. 



of convection zone properties, therefore in terms of IGW char- 
acteristics. The position of the selected points along the evolu- 
tionary track in the HR diagram for the 3 M Q Pop I star is shown 
in Fig.Q] 

Figure|2]shows the temperature at the boundary of both cen- 
tral and external convection zones at the selected evolutionary 
points, which is a way to characterize the depth and extension 
of these regions as the surface and central temperatures (also 
shown in the figure) of the star evolve. In such an intermediate- 
mass star, both central hydrogen- and helium-burning occur in a 
convective core. Wave excitation by core convection on the main 
sequence is discussed elsewhere (Pantillon et al. 2007, Paper III). 
Here we only focus on wave excitation by the surface convection 
zone. 

Figures [3] and [4] illustrate the evolution of the main proper- 
ties that are required to understand the behavior of IGWs in the 
framework described in § 3. Wave excitation is stronger when 
the convective scale 6 C is larger. However, when the turn-over 
timescale 17, becomes too large, this efficiency diminishes. As 
we shall see in § 3 and subsequent sections, the combination 
of these two factors produces large differences in the overall 
efficiency of wave generation as the stellar structure evolves. 
Another important property for wave-induced transport is the 
thermal diffusivity Kj at the top of the radiative region. On the 
main sequence for the star we focus on, Kj is so large that all the 
low frequency waves (here, we mean waves with <x < 3.5 //Hz) 
are dissipated as soon as they are formed, hence the absence of 
a filtered flux in the corresponding models (see § 5). 

Finally, damping also depends strongly on the Brunt-Vaisala 
frequency (Eq. 13, whose thermal part is given by 

N T 2 = ^(Vad- V) (1) 
tip 

with 5 = -(dlnp/dlnT) P4 i, and other symbols have their usual 
meaning. From Fig. @] we see that during the PMS, the star's 
contraction leads to an increase in gravity and a reduction of 
the pressure scale height, which produce a large increase in 
Nt 2 when combined. Later, during the giant phase, the expan- 
sion of the outer layers produces a large decrease in Nj 2 out- 
side the hydrogen-burning shell. More important, core contrac- 
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Fig. 3. Evolution of the surface convection zone properties for the 3 M Pop I evolutionary points from TableQ] It shows the thermal 
conductivity Kj (in cm 2 s ) at the convective boundary, the convective flux F c = pv 3 . (in g cm 4 ) taken 0.5 Hp into the convection 
zone, the mass of the convection zone, the spherical harmonic degree corresponding to the convective scale £ c , the convective 
frequency v c (in /jHz), the typical convective timescale t l = aH P /v c , as well as the total energy luminosity (in g cm 2 s~ 3 ) prior to 
any damping in waves (this is for waves with v < 3.5 fiHz and 1 < ( < 60). 



tion produces a major increase inside this shell; in fact, Nj 2 be- 
comes so large there that (low-frequency) waves will be unable 
to cross that boundary. A large angular momentum deposition 
is expected there, provided there is differential rotation above to 
produce a bias in the local wave flux. 



3. Formalism 

The formalism we use to describe IGW properties is extensively 
described elsewhere (Papers I and II, and TC05). Here, we only 
recall the main features of our model and discuss the critical 
physical principles. 



3.1. Wave spectrum 

In terms of angular momentum evolution, the relevant param- 
eter is the filtered angular momentum luminosity slightly be- 
low the convection envelope (hereafter CE). To get that lumi- 
nosity, we first need to obtain the spectrum of excited waves. 
As we did in previous studies, we apply the Goldreich, Murray 
& Kumar (1994) formalism to IGWs to calculate this spectrum. 
The energy flux per unit frequency Te is then 



™->-5/*£[(£ 



+ £(t + 1) 



xexp[-/i 2 ^+l)/2r 2 ] — 



v L 4 



(cotl) 1512 ' 



(2) 



where £. and [{({ + l)] 1/2 & are the radial and horizontal dis- 
placement wave functions, which are normalized to unit energy 
flux just below the convection zone, v c is the convective veloc- 
ity, L — a MLT H P the radial size of an energy bearing turbulent 
eddy, tl » L/v c the characteristic convective time, and h a is the 
radial size of the largest eddy at depth r with characteristic fre- 



quency of to or higher (h w = L min{ 1 , {2coti) 3 ^ 2 }). The radial 
wave number k r is related to the horizontal wave number kf, by 



N 2 



2 -m= 



N- 



- 1 



W+l) 



(3) 



where A^ 2 is the Brunt-Vaisala frequency. In the convection 
zone, the mode is evanescent and the penetration depth varies 
as y/£({ + 1). The momentum flux per unit frequency Tj is then 
related to the kinetic energy flux by 



Tj (m, (, co) 



2m 



-T E (t,to) 



(4) 



(Goldreich & Nicholson 1989; Zahn, Talon & Matias 1997). We 
integrate this quantity horizontally to get an angular momentum 
luminosity 



£j = 4nr 2 Tj 



(5) 



which, in the absence of dissipation, has the advantage of being 
a conserved quantity (Bretherton 1969; Zahn et al. 1997). This 
angular momentum luminosity, integrated over 6v = 0.2 jjHz, 
is illustrated in Fig. [5] for our selected evolutionary points, that 
have low frequency waves with v < 3.5 yuHz. 

Each wave then travels inward and is damped by thermal 
diffusivity and by viscosity. The local momentum luminosity of 
waves is given by 



-CAr) = ^ -Cj e , m (r C z) exp [-r(r, tr, €)] 



(6) 



crj.m 



where 'cz' refers to the base of the convection zone, t corre- 
sponds to the integration of the local damping rate and takes into 
account the mean molecular weight stratification 



r(r, cr, , 



[{({+!)] 



f 



(K T + v v ) 



N- 
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dr 
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r 5 
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Fig. 5. Evolution of the spectrum of angular momentum luminosity integrated over 5v = 0.2 /iHz (in g s~ 3 prior to any damping of 
IGWs on the PMS (upper row) and from the giant branch to the early AGB in the 3 M Pop I star. Models 4 to 1 1 have no waves 
with v < 3.5 //Hz and are thus not shown here. 



(Zahn et al. 1997). In this expression, Nt" is the thermal part of 
the Brunt-Vaisala frequency (see EqQ}, Kj the thermal diffusiv- 
ity, and v,, the (vertical) turbulent viscosity. Finally cr is the local, 
Doppler-shifted frequency 



cr(r) = a) — m [£l(r) — Q cz ] 



(8) 



and u> is the wave frequency in the reference frame of the convec- 
tion zone. Let us mention that only the radial velocity gradients 
are taken into account in this expression for damping. This is 
because angular momentum transport is dominated by the low- 
frequency waves (cr <k AO, which implies that horizontal gradi- 
ents are much smaller than vertical ones (c/Eq.[3]l. 

When meridional circulation, turbulence, and waves are all 
taken into account, the evolution of angular momentum follows 



d r i„n 1 d r a„ i 13 
07T r z or 



pv v r 



,3Q 



(9) 



(Talon & Zahn 1998), where U is the radial meridional circula- 
tion velocity. Horizontal averaging was performed for this equa- 
tion, and meridional circulation was considered only at first or- 
der. In this paper, this equation is used only to calculate the fast 
SLO's dynamics (see § 13. 2b and thus, U is neglected. The com- 
plete equation is used when longer timescales are involved as 
discussed e.g. in Charbonnel & Talon (2005). 

In this work, we neglect any contribution from convective 
overshoot due to the lack of a usable prescription (although work 



is underway to include this effect, Belkacem et al., in prepara- 
tion). We expect this contribution to be more efficient for low- 
frequency, low-degree waves, and as such, it would have an im- 
pact mostly on the timescale for angular momentum extraction 
in the core (see § 13. 3t . 



3.2. Shear layer oscillation (SLO) 

One key feature when looking at the wave-mean flow interac- 
tion is that the dissipation of IGWs produces an increase in the 
local differential rotation, which is caused by the increased dis- 
sipation of waves that travel in the direction of the shear (see 
Eqs.|7]and[8]l. In conjunction with viscosity, this leads to the for- 
mation of an oscillating doubled-peak shear layer that oscillates 
on a short timescale (Gough & Mclntyre 1998; Ringot 1998; 
Kumar, Talon & Zahn 1999). This shear layer oscillation (here- 
after SLO) occurs if the deposition of angular momentum by 
IGWs is large enough when compared with viscosit)0 (Kim & 
MacGregor 2001; see also a discussion in TC05, § 3.1.). 

To calculate the turbulence associated with this oscillation, 
we relied on a standard prescription for shear turbulence away 
from regions with mean molecular weight gradients 



J^crit 



(rdQ/dr) 2 

Nr 2 



(10) 



4 When the wave flux is too small, one simply reaches a stationary 
regime where the angular momentum deposition by IGWs is compen- 
sated by viscous diffusivity. 
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Fig. 6. Differential wave amplitude variation over a depth of 0.05 R* (see text). 



which takes radiative losses into account (Townsend 1958; 
Maeder 1995). This coefficient is time-averaged over a complete 
oscillation cycle (for details, see TC05, § 3.1). 

The first feature we are interested in is whether or not the 
formation of such an SLO is possible at various evolutionary 
stages, since the associated turbulence could produce mixing in a 
small region below the surface convection zone. Figure|6] shows 
the differential wave amplitude variation in the region of SLO 
formation. The plotted quantity corresponds to 

r L . 2 +0.05R» 

>^dif — ^ ' |AAp r0 g m de ~~ AA retrograde | (11) 

where AA corresponds to the local amplitude variation, caused 
by radiative damping (see Eq. 13 for the £ = m wave and for a 
local gradient of 0.001 yuHz/0.05 R*. This choice of differential 
rotation is quite arbitrary, but TC05 (see their § 3.2) show that, as 
long as differential rotation is not too strong (that is, it remains 
below about 1 /^Hz/0.05 R t ), differential damping is linear. The 
value of 0.05 R* is typical of the thickness of the SLO (TC05). 
When ytdif is large, prograde and retrograde waves are deposited 
at different locations in the region of SLO formation and thus 
contribute to the potential generation of an SLO. As the wave 
frequency rises, the Doppler shift is felt less by the star and ytdif 
remains small. For frequencies that are too low, the waves can be 
damped even before the Doppler shift may take place (see e.g. 
Model 1). 

To discuss the formation of the SLO, ^dif should be multi- 
plied by the angular momentum luminosity of waves, which is 
shown in Fig.[5]fbr the GMK model. This was done for example 
for models 13 and 15 in Fig. [7] For these two models, the waves 
that would give rise to the appearance of an SLO are quite simi- 
lar, as the differential amplitude figure shows (Pig.|6j>. However, 



there is an SLO only in model 13, because in model 15, the wave 
angular momentum luminosity corresponding to such waves is 
simply too small. We note that, for any other excitation model, an 
SLO would be obtained for an angular momentum luminosity on 
the order of 10 33 - 10 34 (or more) in the region with ^dif — 0.01 
to 0.1. Therefore, for any excitation model we can rapidly eval- 
uate whether an SLO should develop from our Figj6] 

The existence of such an SLO has been challenged by Rogers 
& Glatzmeier (2006, hereafter RG06). In 2-D numerical simu- 
lations of penetrative convection in the Sun, they find that IGWs 
are efficiently excited by convection, as did the other authors 
who studied this (Hurlburt et al. 1986, 1994; Andersen 1994; 
Nordlund et al. 1996; Kiraga et al. 2003; Dintrans et al. 2005; 
Rogers & Glatzmeier 2005a,b). In these simulations, wave spec- 
tra are much broader than the ones predicted in theoretical mod- 
els based on Reynolds stresses (Kumar, Talon, & Zahn 1999, 
hereafter KTZ99). However, contrary to the expectations of 
other authors (see in particular the discussion by Kiraga et 
al. 2003), RG06 obtain an energy flu:x0that j s sma n er [ n f] le low- 
frequency regime than the one calculated by KTZ99, and that, 
even though one expects the excitation by convective plumes that 
has been ignored by KTZ99 to yield a large contribution to wave 
generation. 



In their simulation, which has a lower Reynolds number than the 
Sun, but a larger convective flux, the total wave flux is rather large, and 
it scales as 

v 

^waves = F T — (12) 

with F T the total flux, v the velocity of turbulent eddies and c s the sound 
velocity. 
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Fig. 8. Net wave amplitude beyond a depth of 0.05 R* (see text). 



In their simulations, RG06 thus have a much weaker wave 
flux in the low-frequency regime than what is required to gen- 
erate an SLO (see the above discussion on Fig. |7]>, although the 
one for the high frequency waves is much larger. It may thus not 
be surprising that RG06 do not find such an SLO in their simula- 
tion. Furthermore, the size of the SLO being smaller (in KTZ99) 
than the difference between the rotation rate at the base of the 
convection zone (at the equator) and the base of the tachocline, 
an SLO would not appear as a shear reversal in their simulation. 
Last, the RG06 simulation is simply not long enough to show 
the SLO. As pointed out by the authors, in their numerical sim- 
ulation there is differential rotation in the tachocline driven by 
convective plumes. But the total duration of the simulation be- 
ing only one year, it would be impossible to see a superposed 
SLO since its timescale is expected to be of about a few years, 
according to the KTZ99 calculations. Furthermore, since RG06 
have a reduced wave flux compared to KTZ99, this timescale 
would be even longer, and thus clearly beyond the possibilities 
of the RG06 simulation. 



3.3. Filtered angular momentum luminosity 

We must also look at the secular effect of IGWs in the deep in- 
terior (see TC05 for details). In the presence of differential ro- 
tation, the dissipation of prograde and retrograde waves in the 
SLO is not symmetric, and this leads to a finite amount of angu- 
lar momentum being deposited in the interior beyond the SLO. 
Figure|8] shows the net amplitude of waves at a depth of 0.05 R„, 
defined by 



«^net — -^retrograde ^prograde 



for a local gradient of 0.00 1 //Hz/0.05 Multiplied by the an- 
gular momentum luminosity of waves, this is the filtered angular 
momentum luminosity Xy fil given in Table Q] for our selected 
evolutionary points. Let us mention that, in fact, the existence 
of an SLO is not even required to obtain this differential damp- 
ing between prograde and retrograde waves, so as long as differ- 
ential rotation exists at the base of the convection zone, waves 
will have a net impact on the rotation rate of the interior. The 
timescale required to delete the imposed differential rotation is 
of the order of 

t = I6Q.I£./\ (14) 

The SLO's dynamics is studied by solving Eq. (0 with small 
timesteps and using the whole wave spectrum, while for the sec- 
ular evolution of the star, one has to instead use the filtered an- 
gular momentum luminosity (see TC05 for details). All the rel- 
evant quantities are given in Table Q] for our selected evolution- 
ary points. Let us also mention here that differential damping is 
required both for the SLO and for the filtered angular momen- 
tum luminosity. Since this relies on the Doppler shift of the fre- 
quency (see Eqs.|7]and[8]l, angular momentum redistribution will 
be dominated by the low-frequency waves, which experience a 
larger Doppler shift, but that is not so low that they will be imme- 
diately damped. Numerical tests indicate that this occurs around 
(cr * 1 //Hz). 

4. Pre-main sequence 

We start with a fully convective contracting star. As it descends 
along the Hayashi track, a radiative core appears (Fig.Q~|i. In mass 
coordinates, the top of this radiative region migrates towards the 
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Table 1. Selected models along the evolutionary track of a 3 M Pop I star. The filtered angular momentum luminosity £.j is taken 
as 0.05 R below the CZ (or at R cz /2 if it is larger) for a differential rotation of SO. = 0.001 yuHz. Also provided are the star's moment 
of inertia of the radiative zone Ir and the resulting timescale r = IrSQ./X.^j for IGWs. The last column gives the time Af between 
consecutive models (with the convention Af, = f,+i - f,). No waves are generated in models 4 to 11. 
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Fig. 9. Shear layer oscillation in the 3 M Pop I model at the evolutionary points where this feature appears. For each evolutionary 
point, successive profiles are separated by Af = 5 yr. 



surface until the star reaches the main sequence. This is accom- 
panied by a growth of the characteristic convective length scale 
at the bottom of the convection zone i c = 2nr cz /aHp (see Fig.|3]l, 
which temporarily favors the generation of waves (see models 
to 2). 



However, with further retraction of the envelope, the thermal 
diffusivity Kj at the convective boundary increases, which favors 
the disappearance of low-frequency and/or large degree waves 
(see Fig. [5] and Eq. |7). This explains the reduction of the total 
energy flux at the end of the PMS (models 3 to 6) and on the 
MS. 
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0.2 0.4 0.6 0, 

m/M, 

Fig. 4. Evolution of the (thermal part of the) Brunt-Vaisala fre- 
quency vs mass for our 3 M Q Pop I star. (Top) The dashed lines 
correspond to PMS models and the continuous lines, to MS mod- 
els. In PMS models, Nt 2 increases with evolution, while it de- 
creases in MS models. The shade of gray becomes lighter for 
more evolved models. The dashed arrow corresponds to the evo- 
lution during the PMS and the continuous one to the evolution 
during the MS. (Bottom) The continuous lines correspond to 
sub-giants and giants, the short-dashed lines to the clump, and 
the long-dashed lines to early-AGB models. In the outer layers, 
Nt 2 diminishes with time while it is the opposite in the core. The 
arrows indicate time evolution. 
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Fig. 10. Turbulent wave-induced diffusion coefficient at the con- 
vection zone boundary on the pre-main sequence in the 3M 
Pop I star 



migrates outwards with evolution, waves are dissipated more ef- 
ficiently in model 2 than in model 1. In model 3, damping is so 
large that the formation of a SLO is no longer possible, and all 
low-frequency waves are dissipated in model 4. The turbulent 
diffusion coefficient associated with this SLO, D (see Eq.fTOland 
the discussion around it), is shown in Fig. [10] as a function of 
reduced radius. While the shear is much stronger in model 1 
than in model 2, the magnitude of D and width of the turbu- 
lent regions are quite similar. This is because the reduced shear 
is compensated for by the larger thermal diffusivity in model 2 
(cf. Eq. [Tol l. In both models, D is relatively high just below the 
convective envelope, and it drops rapidly with depth. It is neg- 
ligible in the region with a temperature (i.e. T > 2.5 x 10 6 K) 
high enough for efficient lithium burning. This temperature is 
attained at reduced radii of r/R ~ 0.43 and 0.41 in models 1 
and 2 respectively. This implies that, for this mass, the transport 
of elements through wave-induced turbulence has no impact on 
the pre-main sequence evolution of the surface abundance of the 
light elements such as lithium, beryllium, or boron (the latest 
two elements burning at an even higher temperature fj. 

The role of IGWs is more uncertain for the rotation profile. 
In the case of rapid rotation (v > 100 km s~'), the initial rotation 
profile is deleted by meridional circulation in a fraction of the 
main-sequence lifetime. In a slow rotator, however, this is not 
the case, and the memory of angular momentum transport by 
IGWs will influence mixing Detailed evolutionary calculations 
of angular momentum in the presence of internal waves should 
thus take this evolutionary phase into account especially in the 
case of slow rotators (Charbonnel & Talon, in preparation). 



Fig. 7. Differential wave angular momentum luminosity varia- 
tion over a depth of 0.05 R* for models 13 and 15. 



During the pre-main sequence, an SLO exists only for a very 
short period of time. The evolution of the SLO for models 1 
and 2 is shown in Fig. [9] where the successive profiles are sep- 
arated by 5 years. The larger amplitude in model 1 compared 
to model 2 can be understood in terms of thermal diffusivity: 
since Kj increases with radius, as the top of the radiative zone 



5. Main sequence and sub-giant branch 

On the main sequence (models 7 to 9) and early on the sub- 
giant branch (models 10 and 11), the surface convection zone of 
the star remains quite shallow and contains no more than ~ 3 x 
10~ 8 M«. In the bottom part of the convection zone, which is the 



6 No LiBeB depletion is expected on the PMS in a classical model 
for a Pop I star of that mass, that is due to the CE base too cool during 
that phase (see Fig.[2j. 

7 The role of the evolution of the rotation profile on main sequence 
mixing is discussed further in Talon, Richard & Michaud (2006), § 3.2. 
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main driver of IGWs, the mean convective flux is too small (~ 

0. 001 L ) to excite IGWs efficiently (see Fig. [3). Furthermore, 
the thermal diffusivity below the convection zone is quite large, 
and perturbations traveling into the radiative zone are rapidly 
damped instead of becoming waves. This is the case in particular 
for all the low-frequency waves that are needed for driving an 
SLC0(for waves with a frequency v < 3.5 //Hz, Tj - 0). There is 
thus no SLO or any secular effect on the rotation profile from this 
exterior convection zone, both on the MS and at the beginning 
of the sub-giant branch in such an intermediate-mass star. 

This result agrees with the conclusions of Paper I, i.e., the 
total momentum luminosity in waves drops dramatically in main 
sequence Pop I stars with initial masses higher than ~ 1 .4 M , 

1. e. stars originating from the left side of the Li dip, compared to 
stars with lower initial mass. 

Angular momentum transport by IGWs becomes important 
as the star moves farther along the sub-giant branch. Indeed 
the retraction of the convective envelope is accompanied by a 
diminution of the thermal diffusivity and an increase in the con- 
vective flux. To properly assess the importance of waves in the 
star's rotation, we should compare the timescale r (Eq. [T4"l i with 
the lifetime of various evolutionary stages; both quantities are 
given in Table [TJ Model 13 is of special interest here because it 
lies at the end of the sub-giant branch (see Fig. [0 and supports 
a particularly large wave flux, due to a unique combination of 
many factors (including the convective time- and length-scales, 
see Pig . [3j ■ This implies that possible differential rotation, which 
could be a relic of the star's main sequence history and sub- 
sequent contraction, will be strongly reduced by IGWs when 
the star leaves the Hertzsprung gap. This will have a profound 
impact on the subsequent evolution. Let us also mention that, 
close to the core, waves could actually create a large differen- 
tial rotation. This was observed in the numerical simulation by 
RG06 and is also visible when one looks at the Fig. 3 of Talon 
et al. (2002), where, close to the core, strong shears appear in 
the early stages of the simulation. Detailed calculations, includ- 
ing meridional circulation, shear turbulence, and wave induced 
angular momentum transport, are required to come to any con- 
clusion other than that the wave will have a strong impact at this 
stage. 

6. Red giant branch 

As the star evolves to the red giant branch, the surface convection 
zone deepens, and the thermal conductivity below the convection 
envelope is reduced with the photon mean free path (models 12 
to 16, see Fig. [3). Wave excitation during this phase is quite sim- 
ilar to what has been observed in PMS stars. 

A SLO appears at the tip of the RGB (Model 16). The corre- 
sponding diffusion coefficient is shown in Fig.QTJ together with 
the abundance profile of the lightest elements that are affected 
by shell hydrogen-burning. These quantities are plotted against 
SM, which is a relative mass coordinate allowing for a blow-up 
of the radiative region above the hydrogen-burning shell (here- 
after HBS) and is defined as 



on D 



Model 21 log X 



6M 



M r - M H BS 



Mi 



BCE 



-Mu 



(15) 



Here 6M is equal to 1 at the base of the convective envelope and 
at the base of the HBS (where X = 10~ 7 ). The abundance 
profiles at that phase have been modified by dilution during 
the first dredge-up episode and are being reconstructed by shell 




- o 



2 



-6 



Fig. 11. Turbulent wave-induced diffusion coefficient at the con- 
vection zone boundary in the 3 M Q Pop I star (left) on the red gi- 
ant branch and (right) on the early-AGB (wide continuous curve, 
left scale). Superimposed is the logarithm of the mass fraction of 
H, 4 He, 12 C, 13 C and 14 N (right scale). The abscissa is the scaled 
mass coordinate 5M which equals at the base of the hydrogen- 
burning shell and 1 at the base of the CE (see text). Also shown 
is the Brunt-Vaisala frequency (wide dashed curve, right scale). 



hydrogen-burning. We see that D drops very rapidly just below 
the CE. Additionally, the star leaves that phase very rapidly. As a 
consequence, wave-induced turbulence is not expected to mod- 
ify the surface abundances on the RGB for such an intermediate- 
mass Pop I star. 

The rotation profile should hardly be influenced by IGWs for 
that star on the RGB, in view of the timescale t (Eq.[T4Ti which 
is now almost two orders of magnitude higher than the lifetime 
at that stage (see Table[T}. 



7. Clump and early-AGB 

Model 16 marks the ignition of central He -burning in the non- 
degenerate core, and model 17 corresponds to the arrival of the 
star in the so-called clump in the HR diagram. At that moment 
the convective envelope strongly retreats in mass, before deep- 
ening again slowly until the central exhaustion of helium, which 
occurs for model 160. Then the star moves to the early-AGB. 
During that period, the total energy flux in waves remains rela- 
tively modest and there is no SLO. 

The situation changes drastically in model 21 however. At 
that point, the convective envelope has reached its maximum 
depth, which stays almost constant until the first thermal pulse 
occurs (at At » 4 x 10 5 yr after model 21). During that second 
dredge-up, the convective envelope has engulfed most (but not 
all) of the more external step of the 14 N profile left behind by the 
CN-cycle, which explains the steep gradients in the abundance 
profiles of that element and of the carbon isotopes just below the 
convective envelope. In this region, D is relatively high, and we 
expect that the associated mixing leads to an additional slight 
variation in the surface abundance of 12 C and 13 C (but not of the 



We recall that thermal damping is proportional to 1 /cr 4 , see Eq.|7] 



9 The variations in the depth of the CE at that phase are hardly seen in 
Fig. [3] because the total variations represented for log (M cz /M„) cover 
nine orders of magnitude. However, the value for log (M cz /M„) in mod- 
els 16, 17, and 21 is -0.854, -0.138, and -0.747, respectively. 
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carbon isotopic ratio since CN is at the equilibrium), and of 14 N 
with respect to standard second dredge-up. 

For the transport of angular momentum, the rapid increase 
in the Brunt-Vaisala frequency in the HBS could produce strong 
shears that, in complete calculations taking rotational mixing 
into account, would be an important source of mixing during 
the early-AGB. This would modify the composition close to the 
core and would strongly modify the subsequent evolution. 

8. Conclusion 

In this paper, we examined when and how internal gravity waves 
generated by surface convection zones are expected to have an 
impact on stellar models that only have shallow surface convec- 
tion zones during the main sequence. We find that several evolu- 
tionary stages can experience angular momentum redistribution 
by IGWs: 

- During the pre-main sequence, just as the star leaves the 
Hayashi track, waves are efficiently excited and filtered. 
While this has no impact on light element abundances and 
only has a limited effect on fast rotators (see §[4]i, it could 
strongly modify the main-sequence rotation profile of slow 
rotators; 

- At the end of the sub-giant branch (i.e., at the boundary of 
the Hayashi track), angular momentum redistribution is once 
more quite efficient and could totally change the angular mo- 
mentum profile that has been established during the main se- 
quence and the sub-giant phase contraction; 

- During the early-AGB, IGW generation again becomes sub- 
stantial: mixing can be expected to occur both at the con- 
vective boundary, due to the presence of an SLO and at 
the hydrogen-burning-shell boundary, due to the enhanced 
damping of waves caused by the increase in the Brunt- 
Vaisala frequency. 

Stellar models including meridional circulation, shear turbu- 
lence, and wave-induced mixing are underway to quantify this 
impact (Charbonnel & Talon, in preparation). 
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